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Equilibrium charge distribution on a thin, straight conductor has been actively discussed in recent 
years, starting with Griffiths and Li [1], who suspected a uniform distribution, and culminating with 
Maxwell's insightful words [6] as reported by Jackson [5]. The present work addresses the electro- 
static problem of a thin, curved, cylindrical conductor, or a conducting /i/amenf, and shows that the 
corresponding linear charge density slowly tends to uniformity as the inverse of the logarithm of a 
characteristic parameter which is the ratio of the diameter to the smaller of the length and minimum 
radius of curvature of the filament. An alternative derivation of this result directly based on energy 
minimization is developed. These results are based on a general asymptotic analysis of the electric 
field components and potential near a charge filament in the limit of vanishing diameter whereby the 
divergent parts of these quantities are explicitly calculated. It is found that the divergent parts of 
the radial and azimuthal electric field components, as well as the electric potential, are determined 
by the local charge density while the axial component is determined by the local dipole density. For 
a straight filament, these results reduce to those for conducting needles discussed in the literature. 
For curved filaments, the configuration of charges and fields is no longer azimuthally symmetric, 
and there is an additional length scale in the problem arising from the finite radius of curvature 
of the filament. Remarkably, the basic uniformity result survives the added complications, which 
include an azimuthal variation in the surface charge density of the filament. As with the variations 
of linear charge density along the filament, the azimuthal variations vanish with the characteristic 
parameter, only more rapidly. These uniformity results yield an asymptotic formula for the capac- 
itance of a curved filament that generalizes Maxwell's original result. The examples of a straight 
filament with uniform and linearly varying charge densities, as well as a circular filament with a 
uniform charge distribution, are treated analytically and found to be in agreement with the results 
of the general analysis. Numerical calculations illustrating the slow convergence of linear charge dis- 
tribution to uniformity for an elliptical filament are presented, and an interactive computer program 
implementing and animating the numerical calculations is provided. 



I. INTRODUCTION 

A paper of Griffiths and Li in this journal [1] has trig- 
gered a lively discussion on the electrostatic equilibrium 
distribution of charge on a thin, straight, cylindrical con- 
ductor. While suspecting that as the diameter-to-lcngth 
ratio of the conductor approaches zero the charge distri- 
bution on it tends to uniformity at all points sufficiently 
far from the ends, they encountered certain puzzling as- 
pects of the problem in their numerical studies. In par- 
ticular, they identified an apparent contradiction in the 
strict limit of zero diameter, posing the question that 
". . .if the charge density were constant, on a truly one- 
dimensional needle, how could the force on an off-center 
point be zero . . . ?" 

Contributions by Good [2] and Andrews [3] followed 
the work of Griffiths and Li, offering further evidence of 
the uniformity conjecture of the latter, as well as alterna- 
tive methods of analyzing the problem. However, it fell 
to Jackson [4] to provide a clear resolution of the equi- 
librium paradox quoted above by carefully considering 
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the limit of vanishing diameter-to-length ratio and the 
crucial role of the slow, inverse-logarithmic approach of 
charge density to uniformity in that limit. 

A subsequent paper by Jackson [5] reported on his dis- 
covery that the problem in question had already been 
considered by none other than James Clerk Maxwell! 
Maxwell, treating "the electrical capacity of a long nar- 
row cylinder" more than 130 years ago [6], had derived 
the inverse-logarithmic approach to uniformity noted by 
Andrews and highlighted by Jackson, and had eloquently 
summarized the behavior of electrostatic charges on a 
thin, straight conductor. Maxwell's pioneering words on 
the subject turned out to be a very fitting last word as 
well. 

The purpose of this contribution is to investigate the 
distribution of charge on a thin, curved conductor, which 
we will often refer to as a filament. The intuition under- 
lying this inquiry is the expectation that a sufficiently 
thin, smoothly curved filament, when viewed up close, 
would appear long and straight, so that the equilibrium 
behavior of charges on the filament should be nearly uni- 
form just as in the case of straight conductors. This is 
indeed what we find and describe in this paper. 

We approach the problem of charge distribution on 
a physical filament in two steps. The first step is the 
treatment of the fields of any sufficiently smooth, one- 
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dimensional distribution of charge, which we will often 
refer to as an ideal filament, with no regard to its conduc- 
tivity properties. Here we consider the behavior of the 
electric field components and the electrostatic potential 
in the limit of vanishing distance from an arbitrary point 
of the filament, and carry out an asymptotic analysis to 
calculate the divergent parts of these quantities explic- 
itly. We find that the divergent parts of the radial and 
azimuthal electric field components, as well as the electric 
potential, are determined by the linear charge density of 
the filament at the point in question, whereas the axial 
component of the electric field is determined by the spa- 
cial derivative of linear charge density, or equivalently, 
the electric dipole density, at the point. The intuition 
underlying these results is discussed throughout, and the 
results are summarized in Eqs. (19) and ('-'(l). 

In the second step of our analysis, we relate the fields 
of the ideal filament to those of a physical one by cast- 
ing the former as the symmetry axis of a curved, cylin- 
drical conductor. Here, we show that the equipotential 
surfaces of the physical filament approach those of the 
ideal one as its diameter vanishes, and use this to deduce 
the inverse-logarithmic approach of charge distribution to 
uniformity, the primary result of this paper, in Eq. (24). 

The rest of this paper is organized as follows. In §IIA 
we set up the geometry of the ideal filament and for- 
mulate the required regularity conditions. In §IIB we 
present the derivation of the electric field components 
near the ideal filament, and in §IIC we deal with the 
electrostatic potential. We then formulate the require- 
ments of relating the fields of the physical filament to 
those of the ideal one in §IIIA. In §IIIB we consider the 
consequences of enforcing those requirements, whence we 
derive the result expressing the approach of charge den- 
sity to uniformity. An alternative derivation of this result 
using energy minimization is developed in §IIIC. Exam- 
ples of straight and circular filaments providing analyti- 
cal illustrations of Eqs. (19) and (20) are given in §IVA. 
In §IVB we present the results of a numerical study of 
an elliptical filament, including an interactive computer 
program for calculating them. Concluding remarks are 
presented in §V. 



II. ELECTRIC FIELD COMPONENTS NEAR 
AN IDEAL CHARGED FILAMENT 

A. Geometry of a Regular Filament 

Consider a charged filament in the limit of zero diam- 
eter, i.e., idealized as an open or closed curve. We will 
assume this curve to be of class C^, or twice continuously 
differentiable, and parameterized by R(s), si < s < S2, 
where s is the arc length parameter [7]. These regular- 
ity conditions guarantee the existence and continuity of 
the principal normal vector and curvature, as well as the 
tangent vector, at every point. Thus f(s) = d'R{s)/ds 
is the unit vector tangent to the curve at point R(s), 
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FIG. 1: The coordinate system used to describe the geometry 
of a curved, cylindrical filament. The unit vectors {f,O,0) 
form the triad of tangent, normal, and binormal unit vectors 
at point Po of the solid curve. The pair (£>, /3) defines the 
transverse plane at point Po, with the dotted circle delimiting 
the cross section of the curved cylinder. The azimuthal angle 
a and the pair (n, a) define a polar coordinate system on the 
transverse plane. 

with ds = (dR • c?R)^/^. Furthermore, the unit principal 
normal is given by 0{s) = k(s) ^df{s)/ds, where k{s) is 
the curvature at point s. These two unit vectors, together 
with the unit binormal vector /3(s) = t(s) x i'(s), consti- 
tute a local orthogonal triad at each point of the curve. 
Note that the regularity conditions stated above guaran- 
tee that k{s) is a continuous function, allowing us to de- 
fine a radius of curvature g{s) = 1/k{s) which is continu- 
ous and has an infimum which we denote by ~ l/^m- 
Note also that the same conditions imply that the curve 
is non-self- intersecting, i.e., that s ^ s' ^ R(s) 7^ R(s')j 
except for closed filaments for which R(s2) = R(si). The 
minimum radius of curvature and the length of the fila- 
ment are the two relevant length scales for points near 
the filament, and we will find it convenient to use the 
scale parameter L = min(£>m,S2 — si) in the following 
analysis. 

Consider a point Po of the filament corresponding to Sq, 
where si < Sq < S2 for closed filaments and si < Sq < S2 
in case of open filaments. As defined above, t(so), i^(so), 
and /3(so) are the tangent, normal, and binormal unit 
vectors at point Pg . We will also consider a nearby point 
P located at R(so) + eLh{a,SQ) on the normal plane of 
the filament at point Pg, where n(so, a) = cos(a)i5(so) 4- 
sm{a)0{sa), < a < 27r, is a unit vector in the normal 
plane and eL is the distance between Po and P. Here 
and below, e is a dimensionless, positive number which 
represents the distance PgP in units of L. Clearly, as 
e ^ 0, P approaches Pg. 

Figure 1 illustrates the geometry of the filament and 
the coordinate system used to describe it. It also suggests 
that as a ranges from to 27r, the line PgP sweeps the 
surface of a normal disk of radius eL centered on point Po- 
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Relative to this disk, n(a, Sg) is the "radial" unit vector 

at point P and a{a,So) = — sin(a)i>(so) + cos(a)/3(so) 
the azimuthal unit vector. These unit vectors, together 
with t(so), constitute an orthogonal triad at point P sim- 
ilar to the three unit vectors of a cylindrical coordinate 
system. Note that as s ranges from si to S2 and a sweeps 
the interval [0, 27r], the set of points R(s) + eLh{a, s) de- 
fine the surface of a curved, circular cylinder of fixed 
diameter 2eL [8]. 



B. Electric Field Components 

The electric field at any point r not on the filament is 
given by 

E(r) = /c /" 'dsA(s)[r-R(s)]/|r-R(.s)|"^ (1) 

J si 

where A(s) is the linear charge density on the fila- 
ment, assumed to be continuously differentiable, and 
k = (47re(,) . In particular, the electric field at point 
P is given by E[R(so) + eLh{a, Sq)]. 

In the following we will be concerned with the axial, 
radial, and azimuthal components of the electric field at 
point Pq, which we respectively define as the limit of 
E,^{P), E,,{P), and E,,{P) as P approaches P,,. We 
then find 

-E'ax(a, So) = lime^oT(so) • E[R(so) + eLn(a, So)], (2) 
-Bra(a, So) = lime^on(Q;, So) • E[R(so) + eLn(a, So)], (3) 
-Baz(a, So) = lime^oQ:(a, So) ■ E[R(so) + eLh{a, So)]. (4) 

Furthermore, we are only interested in the leading, di- 
vergent contributions to these limits, and will organize 
the following calculations accordingly. It should also be 
noted here that the electric field components at point Pg 
depend on the azimuthal angle a as well as on So. This is 
because the e —^ limit is a singular one, which in gen- 
eral renders the limit dependent on the azimuth of the 
approach to point P„. 

Written in detail, Eqs. (2)-(4) appear as 

£^ax(a, So) = /clinie^o / dsA(s)-r(so) 

•[R(so) - R(s)]|eLn(a, So) + R(s„) - R(s)r^ (5) 

i?ra(Q!, So) = fclime_fO / ds A(s){eL + n(a. So) 

J si 

•[R(so) - R(s)]}|eLfi(a, So) + R(so) - R(s)r^ (6) 

/•S2 

-Eaz(a, So) = fclime_,o / ds A(s)q!(q!, So) 

J si 

•[R(so) - R(s)]|eLfi(a, s„) + R(so) - R(s)r^ (7) 

At this point it is convenient to use the regularity con- 
ditions imposed above on the charged filament, namely 



the C"'^ and properties of A(s) and respectively [9], 
to write 

A(s)=A(so) + (s - s„)A'(so) + A(s, So), (8) 

1 2 

R(s)=R(so) + (s - So)t(so) + 2 ^ «:(so)i>(so) 
+R(s,So), (9) 

where A(s, So) and R(s, So) and are "remainder" terms 
which vanish more rapidly than the terms preceding them 
as s ^ So. In the following, we will use the representa- 
tions of Eqs. ((S) and (9) to organize each electric field 
component into divergent (leading) and finite (remain- 
der) parts as e — > 0. Furthermore, we are primarily in- 
terested in the divergent parts and will not explicitly cal- 
culate the finite remainder terms in the following. We 
note in passing that the remainder terms would in gen- 
eral depend on the detailed structure of the filament and 
its charge distribution. 

We will first consider the contributions to the ra- 
dial component of the electric field near the filament, 

Era\a, So) and Eri\a, So), which arise from the first two 
terms in Eq. (S) respectively. Let us first observe that 
L (the smaller of filament length and minimum radius of 
curvature) is the only physically relevant length scale as 
point Pq is approached. Clearly then, the filament should 
appear as a long, straight line of charge as viewed from 
any point R(so) -l- eLn{a, Sg) if e <C f . Consequently, we 
would intuitively expect Era\P) to approach the elec- 
tric field of an infinite, straight line of constant charge 
density A(so) through point Pg in the direction of t(so) 
as P ^ Po- In other words, £;r^a^(P) ^ 2fcA(so)/eP 
as e 0. Similarly, we expect Erl\P) to approach 
the transverse field of a dipole of moment aA'(so)T(so) 
centered at point Pq, where a is a proportionality con- 
stant (independent of e). However, due to symmetry, 
the transverse field of a dipole located at point Pg van- 
ishes at point P. Thus we expect the leading contribu- 
tion to the radial component of the electric field to equal 
i?r?(a. So) = lime^o2fcA(so)/eL. 

To verify the above expectations for iJi.a(a, Sg), we find 
it convenient to start by making a change of variable from 
s to ^ = (s — Sg)/eL. The resulting expression can then 
be written as 

i;(°)(a. So) = lim,^o[fcA(so)/eL] / ' 

X {I + en(a. So) ■ [R(s„) - R(so + eLO]/eL^} 
x|n(a,So) +C[R(so) - R(so + eLO]/eL^r', (10) 

where Ai = [si — Sg]/eL and A2 = [s2 — Sg]/eL. Recalling 
from Eq. (<!) that [R(so + eL^) - R(so)]/eL^ ^ R'(so) + 
ieL^R"(so) = ''^(so) + ^eL£^K{sg)0{sg) as e — > 0, we can 
rewrite the above result as 

£(°)(a,So)=lim,^o[fcA(so)/eL] / d^[l - ^-eLnisg) 

XCOs(a)^2]|-(^^^^)_^-(^^)|-3^ 
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The first term in the above integral is readily found to 

1 /2 

equal 2 by noting that \n{a,So) — if{so)\ = (1 + C^) 
and that A(i 2) {~oo, +00) in the limit. We there- 
fore find that the first term in Eq. (11) is given by 
lun^^Q2kX{so)/eL, confirming our intuitive expectation 
above. 

The second term in Eq. (11) is a correction to the lead- 
ing contribution wc anticipated, and involves the integral 
Iai '^^^^(l + The asymptotic value of this in- 

tegral is found to equal ln(e~^) by noting that Ai and 
A2 are both of the order of as e 0. Using this 
result, we find that the second term in Eq. (II) is given 
by kX{sa)K{sa) cos(a)lim(:^o ln(e). 

Combining the two contributions from Eq. (11), we 
find 

E^°\a,So) ^ lim,^okXiso)[2{eLy^ + Kiso)cos{a) ln(e)]. 

It is worth pointing out here that just as we anticipated 
the leading term in Eq. (12) as corresponding to a long, 
straight line of charge, we could have arrived at the sec- 
ond contribution by noting that, in the second approxi- 
mation, the filament can be considered a circular arc of 
radius 1/k(so), i.e., the radius of curvature of the fila- 
ment at the point in question. The reader can verify this 
by reference to the treatment of the circular filament in 
§IVA. 

Turning to ii^ra^ (a, Sq) next, and using a procedure sim- 
ilar to that followed above for E^^^a, Sq), we find 

EiJ{a, So = fchme^o / £,d£, — 

X {1 + eft(a, So) • [R(s„) - R(s„ + (M)]/(M] 

X |n(a. So) + e[R(so) - R(so + eLOl/eL^p'. (13) 

Using Eqs. (n) and (0), we implement the e — !■ limit in 
the above expression to find 

/+00 
id^\n{a,s,)-if{s„)\-\ 
-00 

(14) 

which vanishes by symmetry as anticipated. We have 
thus verified that i?i.a(a,So) ^ £'ra(a, Sq) + 2kX{sa)/eL 
as e — > 0, where £'ra(Q!, So) is the finite part of the radial 
electric field at point Pq. 

Our next task is to deal with the axial component of 
the electric field, £'ax(aj Sq). For reasons of symmetry, we 
would expect (a, Sa), the contribution arising from 
A(so) in (8), to vanish since it corresponds to the ax- 
ial electric field of a straight filament of uniform charge 
density as the filament is approached. Indeed following 
a procedure parallel to that used above, we find that 
Ei^{a, So) is given by 

'OO 

x|n(a,So) -C'^(so)^^ (15) 



which vanishes by symmetry as surmised. 

The second contribution to the axial component of the 
electric field, the one arising from (s — So)A'(so) in (S), 
corresponds to the longitudinal field of an electric dipole 
of moment aA'(so)To centered at point Pg, as stated ear- 
lier. Indeed following the procedure established above, 
we find 

£;W(«,So)--fcA'(so)lim,^o / \^d^ 

x|fi(a,So) -^f (so)^^ (16) 

which can be evaluated by using the fact that |n(Q;, So) — 
^t(so)| — (1 +^^)^^^ as noted above. The result is 

(a. So) - lim,^o2fcA'(so)ln(e), (17) 

a logarithmically divergent expression. Wc have thus 
found that the axial component of the electric field 
-Eax(a, So) ^ i?ax(a, So) -l- 2A:A'(so) ln(e) as e ^ 0, where 
Ea.x{a,So) is the finite part of the axial electric field at 
point Po- 

The calculation of the azimuthal component i?az(ctj Sq) 
proceeds along similar lines, with E^\a, So) reducing to 

Ei°\a,s,) = lime^ofcA(so)K(so)sin(a)ln(e-i), (18) 

and Eiz\oi, So) vanishing by symmetry. Note that un- 
like the other components, the azimuthal component of 
the electric field arises from the local curvature of the 
filament, and as such depends on the azimuthal angle a. 
We therefore have the result that the azimuthal com- 
ponent of the electric field i^az(Q:, So) — > -Baz(a, So) + 
A:A(so)k(so) sin(a) ln(e~-'^) as e ^ 0, where i?az(a,So) is 
the finite part of the azimuthal electric field at point Po- 
At this point it is convenient to summarize the above 
results for the electric field components near a charged 
filament: 

The electric field at an interior point Pq of a charged, 
regular, ideal filament is given by 

E(a, So) = E(a, So) + A: lim,_o{2A'(so) ln(e)T(so) 

+A(so)[2(eL)~"^ + k(so) cos(q;) ln(e)]n(Q;, So) 

+A(so) sin(a)K(so) ln(e~-^ )q;(q;, Sq)}, (19) 

where E(q!, So) is the finite part of the electric field at 
point Pq and the notation is as established above. 

It is worth repeating here that the divergent parts of 
the electric field in Eq. (19) arise from the local charge 
and dipole densities at point Po . This is of course what 
is expected physically since, given the regularity require- 
ments on the geometry and distribution of charge on the 
filament [9], there is no other source of divergence than 
the proximity of the charges to the field point. More 
importantly, while divergent parts of the radial and az- 
imuthal electric field components are generated by the 
local charge density, the divergent part of the axial com- 
ponent is generated by the local dipole density. It is the 
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latter feature of the fields and sources that will ensure the 
vanishing of the dipole density (hence approach to unifor- 
mity) for a conducting filament in the limit of vanishing 
diameter, as the argument leading to Eq. (24) will show. 
It should also be emphasized here that the results of this 
section apply to an ideal filament of charge regardless of 
its conductivity properties. 

C. Electrostatic Potential 

We conclude this section by deriving the electrostatic 
potential near the filament. The divergent part of this 
potential is expected to correspond to a long, straight 
filament of constant charge density. Indeed applying the 
methods established above, we find the potential at point 
Pq of the filament as 

$(a. So) = Ha, So) + 2fclim,_oA(so) In(e-i), (20) 

where $(a, Sq) is the finite part of the potential at point 
Pq. As expected, the divergent part of this result is just 
the electrostatic potential of a straight filament of charge 
density A(so), and is independent of the curvature of the 
filament at point Pq. As such, it does not contain any 
information on the aziniuthal component of the electric 
field, or the subleading contribution to the radial com- 
ponent, both of which are present in Eq. (19), while the 
axial and the leading radial components are readily deriv- 
able from it, as will be shown below. The potential terms 
corresponding to the missing components, being of the 
order of eln(e), are finite and vanish in the e — s- limit, 
which is the reason they don't appear in Eq. (20). 

To extract the radial and axial electric field compo- 
nents from the potential, we note that the variables So, 
e, and a on which the terms on the right-hand side of 
Eq. (20) depend are in effect the axial, (rescaled) radial, 
and azimuthal coordinates of a local cylindrical coordi- 
nate system. Consequently, we can write 

V$(a,So) = V$(a,So) + 2fclime^o 
X [f(so)a/aso + n(a, s,)L-^d/de]\{s,) In(e-i). (21) 

Note that while <I>(q!, So) is finite as e — ^ 0, its gradient 
is not and in fact includes a logarithmically divergent 
azimuthal part as mentioned above. 

Implementing the derivatives in Eq. (21), we find 

V(f>(Q!, So) = V$(q!, So) - 2A: lime^o 
X [f(so)A'(so) ln(e) + n(a, So)A(s„)(eL)-^]. (22) 

Note that — V<I'(a, So) ^ E(q;, So), the difference being 
the two a-dependent terms in Eq. (19) which are ab- 
sent in Eq. (20), as anticipated earlier. Taking this into 
account, we find that Eqs. (22) and (19) are in agree- 
ment with respect to the leading radial and axial com- 
ponents. Since the approach to a uniform charge density 
for a thin conductor found in Eq. (24) below is forged 



from these leading components, we are assured that the 
electrostatic potential contains sufficient information to 
imply the same result. This is an important point, since 
the derivation of the uniformity result in §IIIC is based on 
energy considerations using the electrostatic potential. 

Anticipating the uniformity of charge distribution on a 
conducting filament in the limit of vanishing diameter, as 
expressed in Eq. (24) below, we can use the asymptotic 
behavior given in Eq. (20) to derive a universal asymp- 
totic formula for the capacitance per unit length of a 
conducting filament. For such a filament, charge per unit 
length tends to a constant, say A, and the electrostatic 
potential to 2fcA ln(e~-'^), as seen from Eq. (20). Therefore 
the capacitance per unit length, c(e), tends to 

c(e) = [2fcln(l/e)]-\ (23) 

where, it may be recalled from §11, 2eL is the diameter 
of the filament. For a straight filament, Eq. (23) reduces 
to the result first derived by Maxwell [6] ; cf. [4] . 

III. DISTRIBUTION OF CHARGE ON THIN 
CONDUCTING FILAMENTS 

A. Physical versus Ideal Filament 

We are now in a position to consider the physical prob- 
lem of a charged, conducting filament in the limit of van- 
ishing diameter. We will define the physical filament by 
reference to the equipotential surfaces of the idealized 
charge distribution treated in the previous section [10]. 
Specifically, given a value $o for the potential of the con- 
ductor (with respect to infinity), we consider the equipo- 
tential surface S[R(-), A(-), $0] consisting of all points r 
such that ^{y) = $0 to be the surface of the conduct- 
ing, physical filament. Here ^(r) is the (exact) electro- 
static potential produced by the idealized filament. The 
uniqueness property of the Dirichlet problem would then 
assure us that the electrostatic field produced by A(s) in 
the exterior of S is precisely the same as would be pro- 
duced by the physical conductor if held at potential $o- 
Furthermore, the surface charge per unit length of the 
physical conductor would be given by A(s) [11]. 

Having defined the physical conductor by reference to 
an assumed charge distribution, we have shifted the bur- 
den of our analysis to reconciling the resulting profile 
with the actual shape of the physical conductor. In par- 
ticular, we might wonder if, for a given R(s), it is possi- 
ble to choose A(s) so that S[R(-), A(-), $0] has a desired 
form, e.g., that of a curved, circular cylinder of fixed di- 
ameter defined in §IIA. Here we have the benefit of the 
insight gained in the study of straight filaments, espe- 
cially in Refs. [1, 3, 4, 12], where we find an affirmative 
answer. In particular, we know that a conductor with an 
ellipsoidal surface corresponds to a uniform charge dis- 
tribution along the symmetry axis [13]. We also know, 
starting with Maxwell [6] , that a thin, cylindrical conduc- 
tor acquires a near-uniform charge distribution except for 
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endpoints where higher charge densities occur. Note that 
the field and charge configurations of a straight filament 
are azimuthally symmetric. 

In the case of curved filaments of finite diameter, how- 
ever, there's no such symmetry and one cannot in general 
expect a one-dimensional charge distribution to properly 
reproduce the azimuthal distribution of the charges on 
the filament's surface. This shortcoming is of course a 
consequence of the mismatch between the dimensionality 
of an ideal filament and that of the surface of a conduc- 
tor. By the same token, however, any difference between 
the two is expected to be small for a thin filament and 
tend to zero as its diameter vanishes. This expectation 
is in fact born out as the analysis below shows, and we 
will be able to determine the linear charge density of 
a sufficiently thin conductor regardless of its azimuthal 
variations. 



B. Boundary Conditions 

The standard electrostatic boundary condition on a 
conductor is the uniformity of the potential throughout 
its interior and surface, a condition which follows from 
the vanishing of the electric field in the interior of a con- 
ductor in electrostatic equilibrium. Our problem then is 
to constrain A(s) such that S[R(-), A(-), $o] turns out to 
be the curved, circular cylinder of fixed diameter defined 
above. Equivalently, we will require the vanishing of the 
tangential components of the electric field on the surface 
of the conductor, starting with the axial component. 

The axial component was calculated in §IIB where we 
found that i?ax(Q!,So) £'ax(Q!, Sq) -|- 2fcA'(so) ln(e) as 
limj^o- The requirement that i?ax(Q!,So) vanish in the 
limit immediately constrains the charge distribution ac- 
cording to 

A'(s„) = -lim,^o£^ax(a, So)/2fcln(e). (24) 

Recalling that R(so) is an arbitrary point of the fila- 
ment (with the exception of endpoints in case of open fil- 
aments), we conclude that, as its diameter tends to zero, 
the linear charge density on the filament tends to uni- 
formity as the inverse of the logarithm of e. This result 
generalizes and closely parallels those found for straight 
filaments, starting with Maxwell's [6]. 

Next, we consider the azimuthal component of the elec- 
tric field on the surface of the filament. At first blush, 
Eq. (19) seems to imply that the vanishing of the az- 
imuthal component requires the eventual disappearance 
of all charges on the filament! In reality, however, what 
is conveyed by Eq. (1!)) is simply that the equipotential 
surface S[R(-), A(-), $o] generated by an ideal filament 
cannot have a circular cross section except at points 
where the curvature of the filament vanishes. This is 
of course exactly what is expected, since the deviation 
from circularity is needed to balance the electric field 
component generated by the local curvature of the fila- 
ment. This sort of behavior is already familiar from the 



study of straight conductors where variations in the di- 
ameter along the conductor are required to balance the 
effect of unequal distances from the ends, and has been 
discussed in detail in Ref. [5]. Needless to say, this cir- 
cumstance does not imply that circular cross-section fil- 
aments are unphysical, only that our one-dimensional fil- 
ament within is not capable of masquerading them. For- 
tunately, it turns out that the deviation from circularity 
is quite small and disappears rapidly as e ^ 0. This is 
the justification for our use of an ideal charge filament 
to describe the distribution of charge on a thin, physical 
filament. 

To establish the above assertion, we will examine the 
shape of the non-circular cross section generated by the 
ideal filament. The unit normal N(a, So) to the equipo- 
tential surface S[R(-), A(-), $o] at point '[l{so)+eLh{a, Sq) 
is defined by the direction of the electric field in Eq. (19), 
i.e., 

N(a, So) oc lime^o [A'(so) ln(e)f(so) + A(so)(eL)"^ 
xn(a. So) + i sin(a)A(so)K(so) ln(e"^)a(a, Sq)] , (25) 

where we have only included the leading contribution in 
each component. Simplifying further, we find 

N(a, So)=lim£^o{n(a, So) + [LA'(so)/A(so)]e ln(e)f(so) 
+ ^K{so)Le\n{e^^) sin(Q;)Q:(a, So)}. (26) 

Let e = 77(a, Soj^'o) describe the equipotential sur- 
face S[R(-), A(-), $o] in terms of our curved cyhndrical 
coordinates (e, a, Sq). Since the gradient vector V[e — 
77(0;, So, $o)] is normal to the surface in question, its di- 
rection must be the same as the unit normal N(a,So). 
Therefore, using Eq. (^(i), we can write 

d\n[r]{a,So,<i>o)]/da= ^k(so)L 

x?7(a. So, $0) ln[7?(a, So, $0)] sin(a), (27) 

which is valid to leading order as lim^g^oo, or equiv- 
alently, as 77(0, So, $0) — > 0. It is clear from Eq. (27), 
and can be readily verified numerically, that the maxi- 
mum fractional variation in the diameter of the equipo- 
tential surface with respect to a is of the order of 
7;(a, So, 'fo) ln[?7(a. So, $0)], hence quite small for a suf- 
ficiently small diameter [14]. 

It is important to emphasize here that the equipo- 
tential surface e = r]{a, So,^o) is normal to the electric 
field by construction, so that the tangential component 
of the electric field vanishes identically at all points of 
this surface. Consequently, the relevant boundary condi- 
tion in matching this equipotential surface to the surface 
of a curved, cylindrical conductor is the circularity of 
the transverse boundary of the former, i.e., the vanish- 
ing of the partial derivative dhi[a{a, Sq, ^o)]/da. This is 
of course what is guaranteed by Eq. (27) in the limit of 
vanishing diameter. 
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One can also examine the azimuthal variations of 
the surface charge density on the equipotential surface. 
The latter is given by the standard formula eoN(a, Sq) ■ 
E(q;, So). Using Eqs. (19) and (27), we find for the surface 
charge density, 

a(a,So,$„)^^(f!l) {l+l 
27ri77(a,So,$o) L 2 
1 2 
x[-K{so)LT]{a,So,<^„)ln[T]{a,So,<S>o)]sm{a)] |, (28) 

where we have kept the first two terms in 77(0;, Sq, <I>o) and 
suppressed the dependence of the linear charge density 
A(so) on $0. 

The quantity of interest here is dhi[a{a,Sa,^o)]/9a, 
the fractional variation of surface charge density with 
respect to the azimuthal angle. Using Eqs. (27) and (28), 
we find 

d\n[a{a, So, $o)]/9a ~ -k{so)L 

X-q{a, So, $0) ln[l/r/(a, s„, $„)] sin(a), (29) 

to leading order. In other words, the fractional azimuthal 
variations of charge density on the equipotential surface 
vanish as 77(0?, Sq, $0) ln[77(a, Sq, $(,)]. Not surprisingly, 
this is the same rate at which the fractional variations 
of the diameter vanish [cf. Eq. (27)], and significantly 
more rapid than the rate at which the linear charge den- 
sity of the filament tends to uniformity. 

The above results allow us to conclude that as the di- 
ameter of the filament decreases, the azimuthal nonuni- 
formities in diameter and surface charge density disap- 
pear relatively rapidly, confirming our earlier statement 
to that effect. Equivalently, we can conclude that in the 
limit of vanishing diameter, the azimuthal component of 
the electric field on the surface of a thin, circular cylinder 
vanishes as required. 

Let us summarize the above results: 

The density of charge along a regular conducting fila- 
ment of fixed diameter tends to uniformity according to 
Eq. (j-j ), except near endpoints, as the ratio of its diam- 
eter to the smaller of its length and minimum radius of 
curvature tends to zero. 

This is the main result of this paper, and as stated 
earlier, generalizes the result of Maxwell and others for 
straight conductors to curved filaments. As such, it con- 
firms the intuition that the charge distribution along 
a sufficiently smooth, charged filament basically follows 
that of a straight one, and similarly converges to uni- 
formity as the diameter of the filament approaches zero. 
This is so despite the fact that there's an azimuthal vari- 
ation of charge density on the surface of the conductor 
at points of non-zero curvature which serves to balance 
the azimuthal component of the electric field generated 
by deviation from axial symmetry. 



C. Charge Distribution Revisited Energetically 

The tendency of charges in equilibrium on thin, con- 
ducting conductors to become uniformly distributed is 
thus seen to be robust and survive the generalization to 
curved filaments. We arrived at this result in the preced- 
ing section by a detailed consideration of the electric field 
components and equilibrium of charges, together with the 
boundary condition that all points of a conducting body 
must be at the same potential. 

At this juncture it is instructive to reconsider the dis- 
tribution of charge along the filament ab initio using en- 
ergy methods. Specifically, we will look for that distri- 
bution of charge on the filament which minimizes the 
electrostatic energy, given that charges are free to move 
within a conducting medium. 

As a preliminary step, we will provide a proof that the 
state of minimum energy for an assembly of charged con- 
ductors is one of uniform electrostatic potential through- 
out each conducting body [15]. Consider a number of 
charged, insulated conductors in electrostatic equilibrium 
[16]. The charge on the ith conductor will be denoted by 
Qi and the region it occupies by TZj. Since the charges 
are confined to the conductors, we may express the elec- 
trostatic energy of the system as 

^^T. I I d\p{vMv), (30) 

where p(r) is the (volume) charge density and the integral 
is over the indicated region. Furthermore, since the total 
charge on each conductor is given and fixed, we have 

/ d\p{v)^Q,. (31) 
J n. 

Rewriting Eq. (.30) by expressing <1> in terms of p, we 
find 

H^ = V ifc /" d'^rd^r'p{r)\r-r'f^p{r'). (32) 

The problem before us then is to minimize W in Eq. (32) 
subject to the constraints expressed in Eq. (31). Follow- 
ing standard procedure, we associate a Lagrange multi- 
plier fii with the ith constraint and seek to minimize the 
quantity 

^. h [ dVdVp(r)|r-rrV(r') 

-J2.l^d[ d^rp{r)-Q,] (33) 

with respect to variations in the multipliers {pi} as well 
as the charge density p{r). 

Setting the first-order variation with respect to p(r) 
equal to zero gives 

k[ d^r'p{r')\r -r'f^ = <P{r) = Pi, r eTZ^. (34) 
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The analogous equations for the multipHers simply re- 
produce the constraint equations of Eq. (31). Note that 
a r ^ TZi for any i, we get = instead of Eq. (34). 

We conclude from Eq. (^M) that electrostatic energy is 
least when each TZi is an equipotential region, i.e., there's 
a uniform potential throughout the body of each conduc- 
tor. This immediately implies the vanishing of the elec- 
tric field as well as the charge density at all but surface 
points of each conductor. At surface points, on the other 
hand, the equipotential condition implies that the electric 
field is normal to the surface of the conductor. We have 
thus arrived at the familiar results of the electrostatics 
of conductors by energy minimization [17]. 

At this point the argument continues as in §IIIB 
and culminates in Eq. (24), which conveys the inverse- 
logarithmic approach of charge density to uniformity. 
We have thus recovered the main result of this paper 
by means of energy considerations. 

IV. EXAMPLES & NUMERICAL RESULTS 

To illustrate the main results of the foregoing analysis, 
we shall present the cases of straight and circular ideal fil- 
aments as examples of Eqs. (19), (20), and (23) in the fol- 
lowing subsection, and a detailed numerical and graphical 
study of an elliptical conducting filament demonstrating 
the approach of charge density to uniformity with van- 
ishing diameter, Eq. (24), in §IVB. 

A. Straight and Circular Filaments 

Our objective here is an analytical calculation of the 
fields of ideal straight and circular filaments for specified 
charge densities, and a comparison of their asymptotic 
behavior in the limit of e — ^ with those of Eqs. (19) 
and (20). As in the derivation of these equations, we will 
present the limiting behavior as an asymptotic, singu- 
lar part that depends on e, plus a remainder part that 
is finite and independent of e. we will also verify the 
asymptotic formula for capacitance, Eq. (23). 

We shall start with the simplest example of a charge 
filament, namely a straight line of uniform charge den- 
sity A extending from z = —L/2 to z ~ +L/2. Then a 
straightforward calculation gives [18] 

[-(i + + [{i + cy + e^] J 

where ^ = z/L. This result can be used to find the 
asymptotic behavior of the potential as e — > 0: 

$(z) = fcA[ln(l - 4zVi^) + 21ime^oln(e-i)], (36) 

for the electrostatic potential at point z of the filament, 
and 

E(a, z) = 8kXz{L^ - ^z\^ r + lime^o2fcA(eL)"^n(a), 

(37) 



for the electric field, valid for points not too close to the 
ends of the filament [19]. One can verify by inspection 
that these results agree with their general counterparts, 
Eqs. (20) and (19), as well as (23), for the case of a 
uniformly charged straight filament. 

The second example is also a straight filament but with 
a uniform dipolc density, i.e., A(z) = gz, where g is a 
constant. Then we find by straightforward integration 

U(i + c) + [(i + c) +e^] J 

+fc,L{ [(i - o' + [{\ + + (38) 

Taking the limit of this equation as e ^ 0, we find 

$(z) = kg [ln(l - iz^/L^) -2 + 2 lime^o ln(e"^)] z, (39) 

for the electrostatic potential at point z of the filament, 
and 

E(a, z) = kg[8z^/{L^ ~ iz^) - ln(l - iz^/L^) + 2]f 
+2A:.glim,^o[ln(e)f 4- z{eLy^h{a)], (40) 

for the electric field, valid for points not too close to the 
ends of the filament. As in the above, inspection shows 
that these results agree with the general ones given in 
Eqs. (20) and (19). 

Our third example is the case of a circular filament 
of radius R and uniform charge density A, centered at 
the origin of a polar coordinate system (r, 9) . By sym- 
metry, the electrostatic potential of such a filament is 
independent of s = R9 and its axial electric field (i.e., 
the component parallel to r; see Fig. 1) is identically 
zero. On the other hand, unlike straight filaments, the 
potential of a circular one has an azimuthal variation and 
a corresponding electric field component, as will be seen 
below. 

Using the geometry depicted in Fig. 1, we find for the 
electrostatic potential of a circular filament 

$(e, a) = 4yfcA{4[l - e cos(a)] + e^y^^^ 

xK[{[A{l - ecos(a)]/{4[l - ecos(a)] + e^})] , (41) 

where K{-) is the complete elliptic integral of the first 
kind defined as 

K{u)= de[\-usm^{e)]~^''^. (42) 

J 

Note that the azimuthal angle a equals zero at points 
r ~ R — e, < < 27r, and it equals tt at points r = 
R + e, <e <2n. 

The limiting behavior of the potential in Eq. (41) 
can be found using the asymptotic formula K[u) —^ 

1/2 

ln[4/(l — w) ' ], u — > 1. Keeping terms up to order 
e ln(e), we find 

$(a) ^ 2kX{[l + i cos(a)e] ln(8/e) - i cos(a)e}. (43) 
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Up to terms of order ln(e), the right-hand side of Eq. (41) 
reduces to 2fcAln(8/e), in agreement with Eq. (20). Note 
also that the capacitance of the circular filament is given 
by 27ri?A/$, which leads to c(e) [2fcln(8/e)]-^ for its 
capacitance per unit length, in agreement with Eq. (23) 
as e — * 0. 

The electric field of the circular filament is found by 
calculating the negative gradient of Eq. (43). The result 
is 

E(a, 0) = (fcA/i?)| [cos(a)[2 - ln(8)] 
+lime^o[cos(a) ln(e) + 2e"^]] n(a, 6*) 
-sin(a)[l - ln(8) + lime^oln(e)]a(a,6l)|, (44) 

which, given the substitutions L ^ R and n 
agrees with Eq. (19). 

It is instructive to find the azimuthal variations of 
surface charge density for the circular filament using 
Eq. (44). Following the procedure given in §IIIB, we 
find 

d\n[a{a,e)]/da = linie^o^ sin(a)e In(e^i), (45) 

in agreement with the general result in Eq. (29). In 
particular, this result confirms our conclusions in §IIIB 
that azimuthal nonuniformitics vanish relatively rapidly 
as e ^ 0. 

As stated earlier, the three cases treated above serve 
to exemplify the general asymptotic results derived in §11 
and III. 



B. Elliptical Filament 

To illustrate the approach of charge density to unifor- 
mity for a curved filament whose axial field (the electric 
field component parallel to r; see Fig. 1) docs not vanish 
identically due to symmetry, we will consider an ellipti- 
cal filament and numerically investigate the behavior of 
the charge density and electrostatic potential as e — > 0. 
We will also calculate the capacitance of the filament as 
a function of its diameter (thickness) in that limit. 

Consider an the elliptical filament of axes (a, 6), a ~ 
2b, centered at the origin of a polar coordinate system 
(r, 0) with its major axis at 6* = & tt. Using the coordi- 
nate system depicted in Fig. 1, we will numerically calcu- 
late the axial component of the electric field for the first 
quadrant, O<0<7r/2, at the azimuthal angle a = 7r/2, 
with e decreasing from 10~* to 10~^^. Recall that the 
filament thickness is given by 2Le, where in this case 
L = {b'^/a) is the minimum radius of curvature for the 
ellipse which obtains at points = & tt, the points of 
extremal distance from the foci. Note that e = 10^^^ 
corresponds to subatomic dimensions for a filament of, 
say, a = 1 m. 

Starting with the general representation of the axial 
field in Eq. (-5), we make the replacements X{s)ds — > 



1.02 r 




S 0.99 



u 

0.98 ^ ' 

0.0 0.5 1.0 1.5 

Polar angle 6 [rad] 

FIG. 2: Approach of the linear charge density to uniformity 
for an elliptical filament. The three graphs show the charge 
density for the first quadrant of the elliptical filament, at e = 
10"'* (short dash, red), e = 10~* (long dash, green), and 
e — lO^^'^ (solid, blue). Note the logarithmically slow rate of 
convergence to uniformity. 



X{0')d0\ R(so) ^ R(6') = acos{0)x+bsm{0)y+e{b^ /a)z, 
R(s) R.(^'') = acos(6'')i + 6sin(6'')y, where we have as- 
sociated Cartesian unit vectors with the polar plane in 
the standard manner. Note that X{0') = X{s)ds/d0\ 
where ds/d0' = [a^sm^{0') + b^cos^i0')]^/^ . 

The three unit vectors constituting the triad of Fig. 1 
correspond here to t{0) as the unit tangent to the ellipse 
in the coimter-clockwise direction, C'{0) as the inward 
unit normal to the ellipse in the polar plane, and (3 = z 
as the unit binormal, which is a fixed vector in this case. 
Note that we have chosen to look at the potential and the 
T-componcnt of the electric field at the azimuthal angle 
a = 7r/2, corresponding to the filament's "ridge line." 
However, the choice of this angle makes little difference 
for a sufficiently thin filament. 

For the charge density A(s), we assume the form 
Ao + AA(s)/ ln(e) suggested by Eq. (24), and approximate 
AA(s) by a three-term Fourier series. Ideally, the charge 
density would be determined by the requirement that the 
potential $(0, a,e) is uniform, in which case we would 
find an e-dependent result, to wit, AA(s, e). We have de- 
termined the three Fourier coefficients by minimizing the 
variations of the potential at e = 10~^^, and expect the 
potential to be reasonably uniform at e = 10^^^, and less 
so at other values of e. 

The primary objective of the numerical investigation 
is to show that as the diameter of the filament decreases, 
the charge density tends to uniformity while the poten- 
tial remains essentially uniform. The results shown here 
confirm this expectation as the relative variations of the 
charge density and potential are found to be of the order 
of 10~^ and 10"'' for e = lO"'^^. It is worth recalhng 
here that the approach of charge density to uniformity is 
inverse-logarithmic, hence extremely slow. 

Figure 2 is a plot of charge density at three values of e. 
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FIG. 3: Approach of the electrostatic potential to uniformity 
for an elliptical filament. The three graphs show the electro- 
static potential, normalized to unity at ^ = 7r/4, for the first 
quadrant of the elliptical filament, at e = 10^* (short dash, 
red), e — 10~* (long dash, green), and e = 10~^^ (solid, blue). 



namely 10~^ (short dash), 10~^ (long dash), and lO"'^^ 
(solid), for which the relative deviations of charge density 
from uniformity are expected to be in the ratio of 1/4 : 
1/8 : 1/12 respectively. The slow rate of convergence to 
uniformity well known from studies of straight filaments 
is clearly in evidence here. 

The behavior of the electrostatic potential for the same 
values of e as above are shown in Fig. 3, where we have 
normalized the potential to unity at 6 = tt/4. This is of 
course the quantity that must be uniform for a conduct- 
ing filament. Indeed the potential varies about 0.01% at 
e ~ 10^^^ and 0.1% at e = 10^**. Given the approximate 
nature of our fit to AA(s, e), wc find the results shown in 
Fig. 3 in agreement with our expectations. 

It is worth recalling that a strictly uniform charge den- 
sity will in general produce a nonuniform potential along 
the filament regardless of how small its diameter is. This 
nonuniformity is reflected in the axial part of E(q;, So) in 
Eq. (19) for the general case, and exemplified by the first 
term of Eq. {'S7) for the straight filament, both of which 
represent the axial electric field in the limit of e = for 
the case of uniform charge density. That such a resid- 
ual axial electric field persists in the limit of vanishing 
diameter is of course the puzzle raised in Ref. [1] and ex- 
plained in detail in Ref. [4]. Whether straight or curved, 
the residual axial field is caused by asymmetric contribu- 
tions from charges not in the immediate neighborhood of 
a given point and requires an appropriate amount of lo- 
cal dipolc density, AA(s,e), to counteract it. Points that 
receive symmetric contributions from all charges, such as 
the midpoint of a straight filament and points 9 = and 
6 = tt/2 oi the elliptical filament (as well as at the equiv- 
alent points 9 = TT and 9 = 3tt/2 not shown in the fig- 
ures), are therefore exceptions to the general statement. 
As seen in Fig. 3, the potential is stationary at these 
exceptional points for the three values of e represented. 
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FIG. 4; Capacitance per unit length of an elliptical filament 
relative to a circular one as its diameter shrinks by 10 orders 
of magnitude; see text. 



Finally, remembering that Maxwell's original contribu- 
tion focused on the capacitance of the straight filament, 
we show the capacitance per unit length of the ellipti- 
cal filament (normalized to that of a circular one whose 
radius matches the minimum radius of curvature of the 
elliptical filament) in Fig. 4 for e ranging from 10~^ to 
10-^^. Recall from §1IC, Eq. (23), that as e ^ 0, the ca- 
pacitance per unit length of a regular filament approaches 
a universal form that only depends on its diameter. This 
result is born out by the behavior seen in Fig. 4, which 
also echoes the slow rate of convergence to uniformity 
characteristic of the problem. 

We have posted an interactive version of the computer 
program used for the above calculations for use by inter- 
ested readers [21]. It is a code written in Mathematica 
and annotated by explanatory comments. The user can 
manipulate the results or change the parameters, sub- 
ject to the limitations resulting from the nearly singular 
nature of the integrations. 



V. CONCLUDING REMARKS 

The exterior field of an infinitely long, conducting, 
charged, circular cylinder does not depend on its diame- 
ter, and is therefore the same as if its charges were to be 
concentrated on its axis. When finite length effects for 
a thin conductor are taken into account, the correction 
terms to the uniform charge density of the infinite case 
turn out small, except near the cndpoints, and vanish as 
the inverse of the logarithm of the ratio of the diameter 
to length of the conductor [4, 5]. What wc have found 
here is that the same conclusion holds for the charge den- 
sity along a sufficiently smooth conducting filament even 
if curved, albeit with appropriate modifications arising 
from the nonzero curvature of the filament. The intu- 
ition behind this result is the observation that at a point 
where the diameter of the filament is sufficiently small 
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compared to either its length or radius of curvature, the 
filament is locally indistinguishable from a straight one. 
In both curved and straight cases, the leading behavior 
arises from local charge and dipole densities while the 
correction terms originate from charges farther away. 

Beyond the main uniformity result, the two cases dif- 
fer in important details. For example, while the charge 
distribution on any straight filament can be described by 
an appropriate scaling of a universal function, such is not 
the case with curved filaments which can assume an enor- 
mous variety of forms. In particular, a curved filament 
can be closed, in which case complications introduced by 



endpoints would not be present. 

The regularity conditions imposed on the geometry of 
the filament in §IIA are not mere mathematical niceties 
since the uniformity of charge distribution on the fila- 
ment will break down in the presence of a corner or a 
kink on the filament [20]. However, although we have 
excluded such points altogether, it should be noted that 
the uniformity results of this paper are still valid at any 
point of the filament whose neighborhood conforms to 
the smoothness conditions, even if they are violated else- 
where. 
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[19] Specifically, we require that 1/2 — \z/L\ ^ e. 

[20] Mathematically, the controlling scale parameter L, the 
smaller of the length and minimum radius of curvature 
of the filament, would cease to exist, or perhaps vanish, 
at such points, in violation of our regularity conditions 
in §11 A. Physically, a "corner" essentially amounts to a 
point where the radius of curvature of the filament is 
smaller than, or comparable to, its diameter, so that e > 
1, far from the limit required for uniformity. 

[21] The Mathematica 7 program is posted to the URL 
http: //www . csus . edu/ indiv/p/partovimh/pg .nb. 



